A NOTE ON ^-CONVERGENCE 



BJORN GUSTAFSSON AND JACQUELINE MOSSINO 

Abstract. We give a criterion for if-convergence of conductivity ma- 
trices in terms of ordinary weak convergence of the factors in certain 
quotient representations of the matrices. 



1. Introduction 

Questions of homogenization of rapidly varying coefficients in elliptic par- 
tial differential equations have been considered by mathematicians at least 
since the 1970s, and by physicists and engineers much longer (cf. the refer- 
ences in |12j . pQ, A typical example is when the conductivity matrix 
A £ = A £ {x) in an equation 

-div (A £ Vu £ ) = f (1.1) 

in some (bounded) domain C M. n oscillates rapidly at a length scale e > 
and one wants to identify a limiting matrix A = A{x) (presumably less 
oscillating than the A £ ) such that, as e — > 0, the solutions u e converge 
in some weak sense to the solution u of the corresponding homogenized 
equation: 

-div (AVu) = f. 

In the 1970s and 80s, F. Murat, L. Tartar identified the appropriate type 
of convergence, .ff -convergence, for the above type of problems and started 
developing general theories for it. Earlier work was much concerned with 
special cases, like strictly periodic structures. In this little note we give an 
equivalent condition for general ^-convergence in terms of ordinary weak 
convergence for the factors in certain quotient representations of the con- 
ductivity matrices, namely for matrices M £ and P £ appearing when writ- 
ing M £ A £ = P £ . In the special case of stratified media (A £ depending on 
only one of the coordinates) and certain generalizations thereof, explicit 
decompositions of this type have been constructed and used for proving 
.ff-convergence in a series of papers [Sj, [2], |S] The purpose of this 
note is to point out that the existence of such quotient representations is a 
completely general fact in connection with ^-convergence A £ — ► A. 
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The proof consists of an adaptation of methods developed by F. Murat 
and L. Tartar, e.g., in jSj, JH], jEj. ^2], ^01 • In fact, even the result 
can be said to be implicit in their work, but perhaps not explicit. 

For an extension of the results in this note to //-convergence in linear 
elasticity, see 0. 

2. //-CONVERGENCE 

Definition 2.1. Let < a < j3 < oo, f2 C M. n (a bounded domain). Then 
M(a, (3; f2) denotes the set of invertible real-valued nxn matrices A = A(x) 
with entries in L°°(f2) and satisfying almost everywhere in f2 the estimates 

>«iei 2 , 

for £ G W 1 . 

Above the bracket (•, •) denotes the scalar product in M n : (£, rj) = Y^j £,jVj- 
Divergence and curl of vector fields are defined as usual: div D is the scalar 

3 J 

and curl/? is the antisymmetric tensor with components 

(curl-E)i 

We shall also need to take div and curl of matrices, and then the above 
definitions apply to the row vectors, i.e., to the last index. Thus, with 

M = (Mij), P = (Pij), 



dEi dEj 

<).r i <).r; 



(curlM)j jfc 



dMij 8M ik 



dxk dxj 

The parameter e > to be used from now on is by convention restricted 
to take values only in a sequence tending to zero (e.g., e G {1, |, . • • }). 

Definition 2.2. (Tartar U ) Let A, A e G M(a,f3;Q) for some < a < 
[3 < oo and all e > 0. Then A e is said to //-converge to A, 

A e -Z+A, 

as e — > if the following holds. Whenever vector fields -D e , D, E £ , /? G 
L 2 (ft) n satisfy 

£> £ = ^E 5 , (2.1) 

D £ ->> D weakly in L 2 (ft) n , (2.2) 

£ £ E weakly in L 2 (0) n (2.3) 
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with 

{div D 6 } £> o relatively compact in i7 _1 (f2), (2-4) 

{curl E £ } £>0 relatively compact in i?" 1 (^)™ xn (2.5) 

then 

D = AE. (2.6) 

It is well-known 8^ that the iif-limit is unique, that the set M(a, (3; Q) 
is sequentially compact for //-convergence, and also that //-convergence is 
stable under transposition of the matrices: if A 6 — ► A then l A £ — > l A. 



3. The result 

Theorem 3.1. Let A,A £ G M(a,/3;U) for some < a < (3 < oo. Then 

A £ — > A as e — > if and only if there exist n x n matrices M £ , M , P £ , P 
with entries in L 2 (S7) and with M (and hence P) invertible, such that 

M e A e = p E) (3 ^ 

MA = P, (3.2) 

M £ weakly in L 2 (Q) nxn , (3.3) 

P £ -± P weakly in L 2 (fl) nxn . (3.4) 

with 

{curlM e } £>0 relatively compact in H- 1 (ty nxnxn , (3.5) 



{divP £ } £> o relatively compact in H (fi) n . (3-6) 

When this is the case M can be chosen to be the identity matrix I and M £ 
so that curlikF = 0. 

Proof. The proof is based on the "div-curl lemma" of compensated com- 
pactness [5], |14j . We recall that this lemma in general says that if 
f £ ,g £ ,f,g £ L 2 (J7) n are vector fields such that f £ — 1 /, g £ — g weakly in 
L 2 (S7) n and such that div/ £ and the components of curlg £ are all contained 
in a compact subset of H^ 1 ^), then (f £ ,g £ ) (f,g) weakly as distribu- 
tions. 

First we prove the "if" -part of the theorem, which is very easy. So assume 
we have the decompositions (|3.1|) and (|3.2|) with weak convergences M £ — 1 
M and P £ — P as in the statement. We consider E £ and D £ = A £ E £ 
satisfying ((221), J23J), (EH), dSSJ- Then actin § on ^ £ S ives 
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or, in components, 

j j 

Here the div-curl lemma applies for each i and it follows that each of the 
members converge in the sense of distribution, to MD and PE respectively. 
Thus we get 

MD = PE, (3.7) 

which, since M is invertible, is the same as (|2.6|) . 

Now we prove the "only if" part. First we have to construct the matrices 
M £ , M, P £ , P. We may take M = I, P = A. Let A denote the transpose of 
A and let be the i:th unit column vector. Thus e, = Vuj, where Ui is the 
i:th coordinate function: 

Ui(x) — 

Setting also fi = div (*Aej) the equation 

div (UVu) = ft 

is trivially solved by u = u%. 

Now, with /j, m as above (1 < i < n), there is for each e > a unique 
solution it? of the elliptic boundary value problem 



(3.8) 



\ div (A £ Vuf) = f h 
\uf-Ui£Hi(n). 

Using it we define M £ to be the matrix whose i:th row is Vuf. In other 
words, M £ is the matrix with entries 

M £ - = 
« flay 

Then we take P £ to be 

so that t P £ = l A £ (S7u\ , . . . , V« £ ). The so defined matrices M £ , P £ satisfy 

JcurlM e = 0, 
\divP £ = /, 

where / is the vector with components fi € -ff~ 1 (f2). In particular, the 
components of curlM e and divP £ stay within a compact subset of i? _1 (f2). 
From ()3.8|) we get for each 1 < i < n the elliptic estimates 

\\ut\\m(n) < C < oo, 

||U £ VU? \\ L 2 {n) n < C < OO. 

Thus for some subsequence of {e} and some limit fields Vi and <Tj we have 
convergences 

u £ -^Vi, weakly in H l {Vt), (3.9) 
Vuf -» Vvi, weakly in L 2 {Q) n , (3.10) 
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U E Vuf <Ti, weakly in L 2 (VL) n . (3.11) 
The latter convergence together with (|3.8|) shows that 

divo-i = fi. 

At this point we use the mentioned fact that //-convergence carries over 

to the transposed matricies. Thus l A e *A, and since curlVuf = and 
div (M £ ) = fi are compact in // _1 (f2) it follows from the definition of this 
//-convergence that 

o-j =* AVvi. 

Therefore Vi solves the boundary value problem 

div(*AVu i ) = / i) 

But this problem has the unique solution Uj. Thus we conclude that 
Vi = Ui and that =* j4Viij. It also follows that in (|3.9I) - H3.11[) we have 
convergence for the full sequence e (because otherwise one could extract a 
subsequence producing a different solution of (|3.12j0 . With this in mind, 
the convergences (j3.1Uj) . (j3.11j) state exactly that 

M £ — M weakly in L 2 (ft) nxn , 

P £ P weakly in L 2 (Q) nxn . 
This proves the theorem. 

□ 

Remark 3.2. It is clear that the matrices M £ , M, P £ , P appearing in the 
decompositions (j3.1|) . H3.2|) are far from being uniquely determined by A £ , 
A, even when all the conditions (|3,3|) - (|3,6|) are satisfied. For example, none 
of the conditions (|3,l|) - ()3.6j) are affected if M £ , M, P £ , P are multiplied 
from the left by one and the same invertible matrix R = R{x) with bounded 
Lipschitz coefficients, 

For the above reasons one cannot formulate the theorem as saying that 

if (EU, (E2J, (H3J), (E3) hold with M invertible, then A £ A if and 
only if ()3.3|) . ()3.4() hold. However, such a statement is true if an appropriate 
normalization is imposed. Examples of such normalizations are that M = 
I or that P = I. One could also move one of the conclusions to be an 
assumption instead. For example, the following statement is correct (and 
easy to deduce from the theorem): if (|3~T|) . (|3~2|) . (|33|) . (|3~T)|) . (|3~3|) hold, 

then A e — > A if and only if (|3.4|) holds. 

Remark 3.3. From the formulation of the theorem one can easily pass to 
construction of "correctors" (cf. ^3], JH]). Indeed, in order to construct 
correctors for A £ one applies the theorem to the transposed matrices: if 

A £ — > A then l A £ — ► *A. Thus there are matrices l N £ , W, *Q £ , %) such 
that l N £ l A £ =< Q £ , l N l A =*Q, W £ -J N, l Q £ with curlW e and div*Q £ 



relatively compact in H" 1 ^). Here we choose the normalization *iV = /, 
l Q = A. If now we have vector fields D £ , E £ as in the definition of the 
//-convergence (for A £ ) and if an additional weak condition is satisfied, e.g., 
that N £ , Q £ are bounded in L°°(f2) nxn , then the assertion is that 

E e - N £ E -» strongly in L? oc (n) n , 

D £ - Q £ E -> strongly in tf oc (£l) n . 

This means that N £ E and (^-E are good approximations (correctors) of E £ 
and D £ respectively. 

To prove the assertion, first notice that D £ — Q £ E = A £ {E £ — N £ E) and 
that therefore, for every u C f2, the L 2 (u) n -norms of both E £ — N £ E and 
D £ — Q £ E can be estimated from above and below by f (D £ — Q £ E, E £ — 
N £ E) dx. But the div-curl lemma gives that 

(D £ - Q £ E, E £ - N £ E) = (D £ , E £ ) - (D £ , N £ E) - (Q £ E, E £ ) + (Q £ E, N £ E) 

(D, E) - (D, NE) - (QE, E) + (QE, NE) = [D — QE, E - NE) = 

in the sense of distributions (see [7j for further details). From this the 
assertion follows. 

4. Example 

The criterion of //-convergence in the theorem above is particularly useful 
in cases where it is possible to find the matrices M £ and P £ a priori (without 
solving any Dirichlet problem, e.g.). The main example for which this occurs 
is the case of stratified media, i.e., when A £ depends on only one of the 
coordinates, say x±: 

A £ = A £ ( Xl ). 

Then the classical philosophy jSj, JH] is that one should write the relation 
D £ = A £ E £ in such a way that the "bad" components of D £ and E £ are 
expressed in terms of the "good" ones. The good components are those for 
which one has control over the oscillations via the differential equation 
or via compactness assumptions (|2.4|) . (|2.5|) . 

In the stratified case, Df and £fi • ■ • j are g°°d and the rest are bad 
(see ^H] for explanations) , and one thus writes (|2.1|) as (suppressing e for a 
moment) 

D i = ^ i D 1 + J2^ 2 (A ij -^)E j (i>2). 

In order to write this as MD = PE one should take care to multiply the bad 
quantities E\ and D2 ■ ■ ■ D n only by good coefficients, for example constants. 
The simplest and most natural choice of matrices M and P then is 
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/ 1 




P = 




where i > 2 is the row index, j > 2 the column index and Sij the Kronecker 
delta. It is immediate that curlM = 0, divP = 0, and, restoring e again, we 
have that A 6 A if and only if M £ -± M, P 6 -± P weakly in L 2 (Q) nxn . 
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